(N 
> 



Green functions of the Dirac equation with magnetic-solenoid field 

S.P. Gavrilov* D.M. Gitmanj and A. A. Smii-nov* 
February 7, 2008 

-r 

O 

o 

(N 

■ Abstract 

^ , Various Green functions of the Dirac equation with a magnetic-solenoid field 

| (the superposition of the Aharonov-Bohm field and a collinear uniform magnetic 

£Sj ■ field) are constructed and studied. The problem is considered in 2 + 1 and 3 + 1 

dimensions for the natural extension of the Dirac operator (the extension obtained 
from the solenoid regularization) . Representations of the Green functions as proper 
f-^. ■ time integrals are derived. The nonrelativistic limit is considered. For the sake of 

completeness the Green functions of the Klein-Gordon particles are constructed as 
§ i well. 

m 

^ ; 1 Introduction 

In the present article we continue our previous study PJ El Ej of the Dirac equation 
with a magnetic-solenoid field, constructing and studying various Green functions of this 
equation. We recall that the magnetic-solenoid field is the collinear superposition of the 
constant uniform magnetic field and the Aharonov-Bohm (AB) field. The AB field is a 
field of an infinitely long and infinitesimally thin solenoid. Recently the interest in such a 
field configuration has been renewed in connection with planar physics problems, quantum 
Hall effect, and the Aharonov-Bohm effect in cyclotron and synchrotron radiations jHEl 

ei ei mini 

In principle, the Green functions can be constructed whenever complete sets of solu- 
tions of the Dirac equations are available. In this connection, one ought to recall that 
solutions of the Dirac equation with the magnetic-solenoid field in 2 + 1 and 3 + 1 dimen- 
sions were obtained in [Tj. The singularity of the AB field demands a special attention 
to the correct definition of the Dirac operator. The need for self-adjoint extensions in the 
case of the Dirac Hamiltonian with the pure AB field in 2 + 1 dimensions was recognized 
in ^01 \H\ where certain boundary conditions at the origin were established. The regu- 
larized case and peculiarities of the behavior of a spinning particle in the presence of the 
magnetic string were considered in |12| I13j. The problem of the self-adjoint extension of 
the Dirac operator with the magnetic-solenoid field was studied in [21 El • In 2 + 1 
dimensions, a one-parametric family of self-adjoint Dirac Hamiltonians specified by the 
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corresponding boundary conditions at the AB solenoid was constructed, and the spectrum 
and eigenfunctions for each value of the extension parameter were found. In 3 + 1 dimen- 
sions, a two-parametric family of the self-adjoint Dirac Hamiltonians was constructed on 
the condition that the spin polarization is conserved. The corresponding spectrum and 
eigenfunctions for each value of the extension parameters were found as well. In j2l E| the 
procedure of solenoid regularization was also considered. The procedure implies consid- 
ering the finite solenoid and then making its radius go to zero. This procedure specifies 
some particular boundary conditions. The values of the extension parameters correspond- 
ing to the solenoid regularization case were determined in 2 + 1 and 3 + 1 dimensions. 
Further, we call the corresponding extension the natural extension. Nonrelativistic prop- 
agators for the spinless and spin-1/2 particle moving in the pure AB field were considered 
mainly in the relation to the AB effect. The propagator of the spinless particle was found 
in [TJ)1 Uni HZ] as a sum of partial propagators corresponding to homotopically different 
paths in the covering space of the physical background. The nonrelativistic propagator 
of the spin-1/2 particle in the AB field for a particular value of the self-adjoint extension 
parameter was discussed in ^H]- The relativistic scalar case for the AB field was studied 
in PHI- The propagators and the AB effect in general gauge theories were considered in 
[213 121] • Recently, vacuum polarization effects in the AB field have aroused great interest, 
see, for example, |2*2*1 |2*3~] and references therein. 

In the present article, we construct and study the Green functions of the Dirac particle 
in the magnetic-solenoid field in 2 + 1 and 3 + 1 dimensions. The physical importance of 
the problem is stressed by the fact that the knowledge of the Green functions in such a 
configuration allows one to study quantum (and quantum field) effects in the magnetic- 
solenoid field on a regular base. A technical specificity of the problem is related to the 
necessity to take into account all the peculiarities related to the self-adjoint extension 
problem of the Dirac operator in the background under consideration. In Sec. 2 we 
consider the (2 + l)-dimensional case in detail. Here, constructing the Green functions, 
we use the exact solutions of the Dirac equation that are related to the specific values of the 
extension parameter. These values correspond to the natural extension, see above. The 
representations of the Green functions as proper time integrals are derived. In addition, 
we calculate the nonrelativistic Green functions as well. In Sec. 3 we extend the results to 
the (3 + l)-dimensional case. In the Appendix, for the sake of completeness, we present 
the Green functions of the relativistic scalar particle. 

We note that the magnetic-solenoid field belongs to such a type of fields that do not 
violate the vacuum stability. For such fields a unique stable vacuum exists, and quantum 
field definitions of Green functions following below hold true [21]. In particular, the 
causal propagator S c (x, x') and the anticausal propagator S° (x, x') are defined by the 
expressions 



where if)(x) is the quantum spinor field in the Furry representation, satisfying the Dirac 
equation with the magnetic-solenoid field, |0) is the vacuum in this representation. The 
symbol of the T-product acts on both sides: it orders the field operators to its right side 
and antiorders them to its left. The functions S c (x,x'), S c (x,x') can be expressed via 




(1) 



(2) 
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the functions S T (x,x'), 



S c (x, x') = 9 (Ax°) S- (x, x') - 9 (-Ax°) S + {x, x') , Ax° = x° - x ' , (3) 
S~ c (x, x') = 9 (-Ax°) S~ (x, x') - 9 (Ax°) S + {x, x') , (4) 

and the latter can be calculated via a complete set ±ip a (x) of solutions of the Dirac 
equation with the magnetic-solenoid field as 

S*(x,x') ±^a(x) A{X') . (5) 

a 

The solutions with the subscript (+) belong to the positive energy spectrum, whereas 
the solutions with the subscript (— ) belong to the negative energy spectrum. Via a all 
possible quantum numbers are denoted. 

The Dirac equation with the magnetic-solenoid field has the form 

{YPu - M) V> (:r) = . (6) 

Here P u = id v — qA v (x), x = (x u ), q is an algebraic charge, for electrons q = — e < 0, 
M is the electron mass, and A v (x) are potentials of the magnetic-solenoid field. In 
(3 + l)-dimensional case v = 0,1,2,3 and 7^ are the corresponding gamma-matrices. In 
(2 + l)-dimensional case u = 0,1,2 and in what follows, we employ the letter T to denote 
the gamma-matrices. We use for these matrices the following representation 

1 =<7,1 = l<7 , 1 = —1(7 , 



where a 1 are the Pauli matrices. In cylindric coordinates (ip,r), x 1 = r cos ip , x 2 = r sin ip, 
the potentials of the magnetic-solenoid field have the form 

A = 0, eAx = [l + fi + A (r)] eA 2 = - [/„ + n + A (r)] — ^ , 

(A 3 = in 3 + 1), A(r) =eBr 2 /2. (7) 

Here B is the magnitude of the uniform magnetic field, and the magnitude B AB of the 
AB field is given by the expression B AB = $5 (x 1 ) 5 (x 2 ), where $ is the AB-solenoid flux, 
(la + n) = $/$o, $o = 2n/e. It is supposed that Iq is integer and < ji <1. 

The functions S T (x, x') obey the Dirac equation ©, whereas the causal and anticausal 
propagators obey the nonhomogeneous Dirac equations: 

(7"P„ - M) S c (x, x') = -8(x - x') , {YP» - M) S~ c (x, x') = S(x - x') . 

We note that the commutation function S(x,x'), the advanced S (x, x') and the re- 
tarded S ret (x,x') Green functions can be expressed in terms of S c (x,x'), S c (x,x r ) as 
follows 

S (x, x') = S~ (x, x') + S + (x, x') = sgn (Ax°) [S c (x, x') - S~ c (x, x')} , (8) 
S adv (x, x') = -9 (-Ax°) S (x, x') , S ret (x, x') = 9 (Ax°) S {x, x') . (9) 
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2 2+1 dimensional case 



2.1 Sets of exact solutions 

First we study the (2 + l)-dimensional case, for which, as known pOE], the Dirac operator 
with the magnetic-solenoid field in 2 + 1 dimensions possesses a one-parameter family of 
self-adjoint extensions. That provides a one-parameter family of boundary conditions at 
the origin. Following E], we denote the extension parameter as G. Generally speaking, 
the AB symmetry is violated for the spinning particle, which is therefore sensible to the 
solenoid flux sign. As was demonstrated in [21 E], the values = ±7r/2 correspond to the 
natural extension, = —tt/2 if the flux is positive and = tt/2 if the flux is negative. 
Below we present a set of solutions ±ip a (x) of (JB1) which we will use for Green function 
construction according to the formulas We consider the problem separately for two 
values of the extension parameter. 

We start with the case = —tt/2. The positive energy spectrum is given by + e and 
the negative energy spectrum is given by _e , 



+£ = -__£ = VM 2 + u . (10) 

Both branches are determined by the spectrum of the quantity to which is defined be- 
low. The solutions ±ip a (x) can be expressed via the solutions u (x) of the squared Dirac 
equation. The latter solutions have the form 

±u m ,i,a (x) = e' l±ex °u mM (x±) , 

x±= (x\x 2 ) , m = 0, 1,... , Z = 0, ±1, ... , a = ±l, (11) 



where 



%,o,+i ( x ±) = V^9o (<p) <Pm,o,+i (r) V +1 , 

u m ,o -l (x±) = ^jg (<p) C-i ( r ) v -i , 7 = e \B\ 



and 



9i(<p) = ^= ex P { i( P 

V Z7T 



/-/ -i(l + a 3 ) 



n / o 



The functions 4> m< i,a- ( r ) , m -i( r ) are expressed via the Laguerre functions I m + a ,m(p) as 

<t>m,lA r ) = tm+Wlm (p) , m ,_].(»") = im-^m (p) , 

p = 7 r 2 /2 , z/ = fi + I - (1 + a) /2 . (12) 

We recall that the Laguerre functions I m + a ,m{p) are related to the Laguerre polynomials 
L m (x) (8.970, 8.972.1 [25]) as 



Jm+Q - m(x) ~ y r (m + a + i) e x/2xa/2L -( x ) 
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For the magnetic field B > 0, the spectrum of uj corresponding to the functions u m ^^ (x±) 
is 

f 2 7 (m + Z + /i), Z-(l + <r)/2>0 , . 

\ 2 7 (m+(l + a)/2), / - (1 + a) /2 < ' 1 j 

except the functions w m ,o,-i (#±) for which the spectrum of is 

cu = 2 7 m. (14) 

Then the complete set ±ip a with a = (m, Z) has the form 

± Vw (*) = iv (rp + M) ±1^,,.! (x) . (15) 

The latter form provides correct expressions both for uj ^ and uj = 0, since the states 
with uj = can only be expressed in terms of the spinors with a = — 1 (we note that 
+ ifj = for u — 0, nevertheless it is convenient to remain in (jllj) m + with + e = M). The 
normalization factor with respect to the usual inner product (t/},ip') = J ift '(x)ip'(x)dx 

[2\±e\(\ ± e\-M)}- 1/2 , u ± , 
[2M]" 1 , cu = 

The quantum number Z characterizes the angular momentum of the particle, m is the 
radial quantum number, see PQ. 

For B < the spectrum of states differs nontrivially from the expressions given by 
Eqs. (fT3j) and (JHJ). Here uj corresponding to u m ^ a {r) is 

j2 7 (m-I + l- / /),I-(l+ ( 7)/2<0 

\ 2 7 (m+(l-«r)/2), Z-(l + <x)/2>0 ' 1 ' 

except the functions u OT) o,-i (;£_l) for which the spectrum of uj is 

cj = 2 7 (m + 1 - fi) . (17) 

Now we go to the case with the extension parameter = n/2. We recall that one needs 
for self-adjoint extensions of the radial Dirac Hamiltonian only in the subspace I = to 
which we refer as to the critical subspace. Thus, the only solutions in the / = subspace 
must be subjected to the one of asymptotic condition from a one-parametric family of 
boundary conditions as r — > 0. By this reason for G = 7r/2, the solutions only differ from 
(fTTJ) in the subspace 1 = 0, 

U m ,0,+l(x±) = ^g (<p)(j) l ^ +1 (r)v +lj C,+l( r ) = J m+/i-l,m (p) , 

u m ,o,-i (x±) = y/ig Q (tp) 4>m,Q,-i (r) u_i , (18) 

where the spectrum for M m ,o,+i (x±) is given as 

uj = 2 7 (m + /i) , B > 0, (19) 
uj = 2 7 m, B < . (20) 



5 



2.2 Construction of Green functions 



The main point in constructing the Green functions is the summations in the representa- 
tion (jHJ). In the case under consideration, this summation can be done with the help of 
special relations which can be established for the solutions of the Dirac equation. 

Let us start with the calculation of the Green functions for the extension parameter 
= — 7r/2 and B > 0. In this case, taking into account that the eigenfunctions u of the 
equation [(TPi) 2 + uA u = corresponding to any w ^ obey the equations 

TP ± ±u m+)l _ a (x) = -iy/uJ ±u m _^ a (x), I < , P± = (0, Pi, P 2 ) , 

TP± ±u m ,J -<r (x) = iVuj±u m> i )a (x) , I > 1, m± = m + (1 ± a) /2 , (21) 

and the explicit form of the solutions ±ip m ,i, one can verify that for \e\^ M the following 
relations hold true, 



±VW (X) ±1p m ,l ( X ') 



2 ± e 

(rp + M) J_ e -*±^- 

2 ± £ 



cr=±l 



er=±l 



Z > 1, (221 



where 



Aip = (p — <p' , a 



JL e <P-'o-(l+<r)/2]A V j f fl W f fl ') 

„ ° - I m+o,m Vr/ ± m+a,m \y j 

Aix 



(23) 



fi + l-(l + a)/2, l>\ 
- \p. + I - (1 + a) /2] , Z < 



L tl = (l±a 3 ) /2 



The above relations and Eqs. (JHJ), © allow us to represent the causal Green function in 
the following form 



S c (x, x') = (TP + M) A c (x, x') 



A c (x,x') =iJ2 



e (Ax ) 



2+e 



e (-As ) 



-i _eAx° 



2_e 



(f>m,l,a {X±, x' ± ) E a .(24) 



Then we can use the representations 



-i-eAx 



i) 



9 (Ax°) 



2+e 



9 (-Ax ) 



2_e 



2m 



e -ip Ax° 



-dpo 



e 2 — pi — it 



e -i^-p 2 ) s ds 



(25) 
(26) 



in Eq. (|24jl . Integrating over p > we obtain finally 



A c (x,x') 



f (x, x , s) ds 



*JLjn/4-iM 2 s 



-pre 4» e ' e "2 

2 M' /2 



l ,i )0 -(x_L,x / j_)H (7 . (27) 



The path of the integration over s is deformed so that it goes slightly below the singular 
points Sk = kir/j, k — 1,2, ... . 
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Using JHJ), (|221)) and the representation 

_ ^ (_Ax°) + 9 (Ax°) = — / -4-^ d Po 

—00 

2 m 2 



instead of (f2l)|). (}2T)j) we obtain from (jlj) 



/— 00 
f{x,x',s)ds, (29) 
-0 



Z7T 

> 



where / (x, x', s) is given by Eq. (|27|) . The negative values for s are defined as s = \s \ e 
and the path of integration over s is deformed so that it goes slightly below the singular 
points — Sfc. 

We now consider the summations in (|27jl. Applying the formula (8.976 (1) |25| ) we 
can sum over m to get 

^ e~ i2mJS I m+a , m (p) /m+a,m (p') = eX P \ % - (p + f/) COt ( 7 s) I 
m=0 ^ ' 

e ia~fs e iis _ i2l2L 

x o . ■ , — re 2 J a (z), z= y/pff/ sin (7s) , (30) 
2z sin (7s) 

where J a (z) are the Bessel functions (8.402 [25]), and for negative s we take args = — n+0. 

Similar results can be obtained for the case B < 0. Here one should use the solutions 
corresponding to the spectrum of uj ()16p . (JT7J). Then these results can be united to obtain 
expressions which hold true for any sign of B, 

00 

/ (x, x', s)=J2h (*> x 'i s ) > fi (*> x '> S )= A (*) E $ ^ ( s ) z~ l ° eBs ^ > 

/=— CO (T=±l 



A (S) = 87T3/V/2 sin (eBs) 6XP iT " ^ " ^ 

x exp J- ' (A:Eo)2 + — (r 2 + r' 2 ) cot (efls) 1 , (31) 
4s 4 

$ i, a (s) = e^e^'+^'e-^^J^+A (z) , l a = I - (1 + a) /2, Z ^ , 

Now we consider the summation over I. One can see that the following relations hold 
true 

00 00 

£ (s) = £ (s) = e-^ Bs Y (z, Aif - eBs, p) , 

1=1 1=1 

—00 —00 

(«) = £ ®*+i>+i ( s ) = e^ eBs Y (z, -Aip + eBs, -fi) , 



Z=-l {=-1 
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where 



Y (z, r], n) = ai (z) + Y (z, ri,fi), Y (z, 77, fx) = ^ 



a, (z) , a x (z) = e ir > 1 (-i 



1=2 



(33) 

The evaluation of the sum in (}3*3*j) can be done in a similar way to what was done in . 
There exist all d z a\ (z) on the half-line, < z < oo, and the relation (8.471 (2) |25j), 
d z J u (z) = [J u -i (z) — J u +i (z)) /2, can be used. The series Y (z,r),fi) converges and the 
series of derivatives Yli=2®z a i ( z ) converges uniformly in (0, oo). It is sufficient condition 
to write down d z Y (z,7),fj,) = Yl*i=2®z a i ( z )- Thus, one arrives to a differential equation 
with respect to Y (z, r], //), 



(34) 



(35) 



j-Y (z, 1], fx) = -Y (z, 7}, fx) i cos 7] + ~ (-if [-ie lv J u (z) + J 1+fM (z)] . 
that is true on the half-line, < z < oo. The solution of (jMj) reads 

Y (z, V,V) = \ Hf f [-ie^J, (y) + J 1+/4 (y)] dy . 

This is also valid for Y (z, —77, —fx). 

It is useful to introduce the following function 

f nc (x, x', s) = 53 /i (x, x', s) . 

It defines the part of the Green functions that is the same for all extensions. With the 
help of the function Y (z, rj, fx) (JSHJ), (|33|) one can write 



f nc (X, X', S)=A (s) e -^Bs e -ieBsa |y ^ ^ _ ^ ^ + y ^ + ^ _ ^ 



+ 



e ^ J^( z ) 



-i(Aif-eBs) 



iTriX — ji) 



2 Jl-u (Z) 



— +1 



(36) 



The function / (x,x', s) is specific for each extension. It is reasonable to mark it with a 
superscript that assumes the values of the extension parameter. Thus, for 6 = —n/2, 



(-7T/2) 





e~ lAip e 



iTr(l-fi) 



2 (^) 



+ e 



-ieBsa 3 2Lf» T / \ 



Accordingly, the function / (x, x', s) acquires the same superscript, 

/ (x, X , s) = fnc (x, X , s) + f ^ (x, X , s) 

For the extension parameter G = tt/2, one obtains 
ft /2) (x,x',s) = A(s)e-^ Bs 



(37) 
(38) 



■ a iir(u,— 1) 

e-^e— ^J M _! z ~ +1 + e" 



(*, x', S ) = f nc (x, x', s) + ff™ ( X| ^ s ) 



(tt/2) 



(39) 



Besides, one can consider particles with "spin down" polarization in 2 + 1 dimensions. 
The corresponding wave functions ij)^~ l > (x) can be presented as 

^G 1 ) ( x ) = a l (TP — M) u (x) , 
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where u (x) are solutions of the squared Dirac equation. The propagator related to 
such particles can be expressed in terms of the function A c (x, x') (j24j) . 

5 ( c _ 1} (x, x') = -a 1 (TP - M) A c (x, x') a 1 . 

At this point we ought to make some remarks. 

One can see that there exists a simple relation between scalar Green functions and 
Green functions of the squared Dirac equation (for the above considered extensions). 
Consider this relation in the example of causal Green functions. First of all, we note 
that the Klein-Gordon equation differs from the squared Dirac equation by the Zeeman 
interaction term. Then we can see (remembering the origin of the quantum number 
I for both spinning and spinless particles) that the scalar propagator can be derived 
from A c (x, x') by only retaining the terms with a = — 1 only. The term eBa 3 , which is 
responsible for the Zeeman interaction with the uniform magnetic field, has to be removed. 
The Zeeman interaction with the solenoid flux, influencing the terms with 1 = 0, depends 
on the flux sign and can be repulsive or attractive. The repulsive contact interaction 
case is physically equivalent to the spinless case, since in both cases the corresponding 
wave functions vanish at the origin. The necessary boundary condition is realized for the 
extension parameter = n/2. Thus, one can obtain the scalar Green functions using 
the coefficients of E_i in fi (x, x', s) (|31|). and (x, x', s) ()39j). By following such 

prescriptions, one arrives at the expression (|5H|) obtained by direct calculation. 

In the spinless case there is no physically preferred orientation of the plane x x x 2 . 
Therefore, the solenoid flux direction does not matter, i.e., the AB symmetry, l — > l + 1, 
is conserved. The direction of the uniform magnetic field does not matter as well. This can 
be observed from the explicit form of the Green functions (|58j) where the change B — > —B 
is equivalent to the choice of the opposite orientation of the plane, I — > —I, Atp — > —Aip, 
$ — > — $. In the spinning case the given spin direction breaks the symmetry related to 
the plane orientation. The Zeeman interaction of the spin with the background violates 
the AB symmetry as well as the symmetry with respect to the change B — > —B. 

As is known, influence of the solenoid flux on the particle is observed only when the 
flux is not equal to an integral number of quanta (/i 7^ 0). In this connection it is 
instructive to consider the Green functions for the particular case \x = 0. We note that 
the part f nc (x, x', s) (J5fijl of the function / (x, x', s) is regular everywhere, while the part 
fo (x,x',s) is singular at the origin. Thus, taking the limit \x — > in and using the 
relation J x (y) = — J' Q (y) we get, 

f nc (x, x', s) = A (s) e" ieW { e — s(A^eB S ) _ Jq {z) + [ Jo (z) + . e -i(A v - eBs)ji s+i J 

The corresponding expression for / ( be obtained in the following way. We 

restrict the range ofzto0<5<z<oc, where 5 <C 1. Then we take the limit /1 — > 
and use the continuity of the Bessel functions with respect to its index. At the end we 
construct the analytic continuation of the obtained expressions over the interval (0,5). 
Thus, starting from either (|S7jl or (J39~j) we get, 

fo (x, x', s) = A (s) [-ie-^J, (z) E +1 + e ieBs J (z) S_ x ] , 
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where the superscript is no longer necessary. Thus, the explicit form of / (x, x', s) is 

„ , . , eB J m % (AxqY . , ,.> . r> , 
!{X ' X • S) = 8^/V/2 Bin(eBs) ™ P |T~ — tM " ~ ,eBs " 

x exp { - fl, A y + *g (r» + r") cot («B.) - ^ ~ eBs > } . (40 ) 

Making a transformation to Cartesian coordinates in ()4Uj) and setting Z = 0, one can 
obtain the known result of the uniform magnetic field, see for example (2ZJ. 

2.3 Nonrelativistic case 

Consideration of the Green functions in the background under question in the nonrelativis- 
tic case is important for various physical applications. Below we study this case in detail. 
The solutions of the Schrodinger equation for "spin up" particles (+) and antiparticles 
(— ) in the case 6 = — n/2 read, 

+ <t> m ,l(x) = e-^^e*< | -''- 1 ^ mA+1 (r) lJ B=^ I (41) 

iEx° / 7 -i(l-l 



-6n,o(*) = e-^°^e*^C,-i(r), (42) 

where the values uj m ^ a are defined by m,l,a with the help of formulas (JUJ), (JHJ) for 
B > 0, and (fTIj|). (|T7jl for _B < 0. The solutions +0 mj z (x) ( _0 mi ; (x) ) for the "spin down" 
case can be obtained from the solutions -4> m ,i (x) ( +4> m ,i (x) ) for the "spin-up" case with 
the change if -> -</? in (gTJ), (gU). 

The retarded Green functions for particles and antiparticles are defined as 

S-*.(±) ( x , x') = 5 (As ) S l ±] x ') . Sf *0 = < E (*) ±<i » 

*') = £<Sf ±) (s.*') > ( 43 ) 

where the part (x, x') is the same for all extensions, whereas Sq^ (x, x') is specific for 
each extension. Carrying out the summations in (J4HJ) one obtains 



Sl ±] (X, X>) = A nr (X, x') e Wr e ± i fe-^)A V , e -^ ±+M | 7 . e -Hi±±H ^ 

(44) 
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A nr (x, x ) = - — — — r exp 
4-7T sin (7T) 



- (p + p') cot ( 7 r) 



^ (*, *') = A nr (X, x') e-^o^d+Mjeflr ^ ( ^) _ e^^-^ 

+y (z nr , A<p - e£r, p) + Y (z nr , -Acp + eBr, -p)} , (45) 

tfH ( Xj X ') = A nr (X, X') e *A Ve i(l-M)eBr 

x {F (z nr , -Aip - eBr, p) + Y (z nr , Aip + e^r, -p) } , (46) 
z nr = ypT 7 / sin (jr) , r = Ax°/2M , Z ± = / - (1 ± 1) /2, / ^ , 
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whereas for 1 = 0, 



A nr (x, x') e-Co+DA^-^^ J±(i _^ ^ ; (4?) 
A nr (x, x') e^e^-^e^J^ (z nr ) . (48) 

The Green function in the " spin down" case can be obtained with the change Aip — > —Af 
in (jlH)-(jiHJ) and with the change 5 (±) by in all the functions S(x,x') in 
Thus, one can see that the Green functions for the nonrelativistic particle is irregular at 
r = when the contact interaction is attractive. 

We note that for the limiting case B = (the uniform magnetic field is absent) 
Sj ^ (x,x') coincide with the known expression for the spinless particle |15 | IT6 | ITT], 
which is natural in the case of a repulsive contact interaction. While Si + ^ n ^ (x, x') for 
B = coincide with the corresponding expressions obtained in the paper |18j . 



o(+)(=pt/2) r T T i\ 

cH(T»/2) , m 



3 3+1 dimensional case 

To obtain the Green functions in 3 + 1 dimensions we use the orthonormalized solutions 
±^ P3 ,m,i,a (x) of the Dirac equation found in [3 El- The quantum numbers m, I have the 
same meaning as in the (2 + l)-dimensional case, ps is the a; 3 -component of the momen- 
tum, and a is the spin quantum number. The positive energy spectrum is given by + e and 
the negative energy spectrum is given by _e . They both are expressed via the quantity 
u as 

+e = --e= ^M 2 +p 2 3 + u. (49) 

The spectra of uj are given in (JJ3J), (fT3)l for B > 0, and in (JTSj). (|T7|) for B < 0. For u; 7^ 0, 

one can present the solutions ±^ P3t m,i,a in the following form, 

±Vp 3 , m ,l,a (X) = N (YPu + M) ±U p , :m ^ a (X) , 

±u P3>mM (x) = -L e -*^°-^ 3 Um ^( x± ) ; 

*W = f {X (r ) ) ' N = I 2 l ±£ l ^ ±£ l + P3 )]" V2 ' ^ 

whereas for c<j = 0, 

± ^ P3i0 ,/,- f (x) = N (>fP v + M) ±U P3 ,o,i,^ = s g n ( 5 ) , 

where u m ,i,a (x±) are the two-spinors defined in (fTTjl . 

We are going to construct the Green functions using the solutions that correspond to 
the natural extensions of the Dirac operator, i.e., for the extension parameters chosen as 
0+i = 0_i = 0, and = ±7r/2. First we consider the case = — it/2, and B > 0. We 
note that for u 7^ 0, 

r r v P± u U rri7 i _ CT = iy/uU m> i ta , I > 1 , 

l u P±. v U m+ ,i -a = -iVuU m _^ a , I < , (51) 
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where P± = (0, Pi,P2,0). The summations in (j5J) can be done similarly to the (2 + 1)- 
dimensional case by the help of some important relations derived by us for the solutions 
((HQ]) . Namely, for the states with a given the following relations hold true 

(7=±1 

= E 2^ ^ + M ) 5 + (tS3 ) ±^.«M,«t (a:, x') , I > 1 , 



<T=±1 



cr=±l 

^ 2^i (7 ' P " + M) ^ t 1 + ^ i^.m+^-a (ar, x') , / < 



=±i 



(52) 



and for uj = 0, we have 

— 1 1 

±*Pa,o,i,-i (&) ±*P8,o,i,-i 0*0 = ^ (V^/ + M) - (1 - £ 3 ) ±^3,0^-1 (x,x') , 

where 



-i ±eAx°— ip3 Ax 3 



± ( Pp3,m,l,cr {X,X ) 27i- e 



0m,i,<7 (x_L, x' ± ) , Ax 3 = X 3 - X 
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The functions 4> m ,i,a (x_i_, defined in (|2Hjl . Therefore, 

S c (x,x') = + A c (x,x') , 

/OO 1 
dp 3 (l +(T £ 3 ) 
. . _ -oo * 



m,L. 



9 (Ax°) - — + P3 , m ,/ iCT (x, a/) - 6> (- Ax°) 7^ -(j) P3l m,i )(T (x, x') 



2_£ 



(53) 



(54) 



Applying the relations ([25)1 . (|26|h one obtains the proper time integral representation for 
A c , ' 

/•oo 00 

A c (x,x') = / f(x,x',s)ds, f(x,x',s)= Y] fl( x , x ', s ) , 
/, (x, x', s) = D (s) 00 ^ i<TeBs ~ (1 + ^S 3 ) , 



=±1 



£>(s1 



167r 2 s sin (eBs) [4s 



exp - [(Ax 3 f - (Ax ) 2 ] -iM 2 s 



x exp <{ -il Aip + — (r 2 + r' 2 ) cot (eBs) }> , 



(55) 



where <&i a (s) are defined in (|32|). 
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Carrying out similar calculations for B < one can verify that (|55|) is valid for both 
signs of B. Therefore, for any sign of B, we get 

f nc (x,x',s) = ^2fi (x,x',s) 

1^0 



D {a) e ~ ieB <^) {Y (z, Ay? - eBs, /i) + Y (z, -A<p + eBs, -/i) 
+ 



e 2 J^{z) 



< A ^e-^J 1 .,(z)] 1 -(l + ^) 



ft 



(-7T/2) 



-D (s) e~ ifieBs 
2 v ' 



e~ iAip e 



i7r(l — fi) 



Ji_„ (z) (1 + S 3 ) 



+e ieBs eTj_,(.) (l-£ 3 ) 
/*■ ^ ^ (a;, x , s) = / nc (x, x , s) + /g (x, x , s) . 
Using the corresponding solutions for the case 6 = ir/2, we obtain 



(56) 



ft /2) (x,x',s) = ^D(s)e-^ Bs 
+e- Bs e-^J,{z){l-Y?) , 

^ (x, x , s) = / nc (x, x , s) + /g ^ (x, x , s) 



1 



i7r{(i — 1) 



J M _! (2) (1 + S 3 ) 



(57) 



4 Summary 

Various Green functions of the Dirac equation with the magnetic-solenoid field are con- 
structed as sums over exact solutions of this equation. We stress that doing that we had to 
take into account all the peculiarities related to the self-adjoint extension problem of the 
Dirac operator in the background under consideration. Both 2 + 1 and 3 + 1 dimensional 
cases are considered. Compact form for the Green functions was obtained thanks to the 
important relations (|2^j) and (fo^j) derived by us for the exact solutions under considera- 
tion. The representations of the Green functions as proper time integrals are constructed. 
The kernels of the proper time integrals are represented both as infinite sums over the 
orbital quantum number / and as simple integrals. The Green functions are obtained for 
two natural self-adjoint extensions, one for the positive solenoid flux and the other one for 
the negative solenoid flux. The physical motivation for the choice of these extensions is 
their correspondence to the presence of the point-like magnetic field at the origin and their 
close relation to the MIT boundary conditions [213 12H1 12H] • Thus, the considered cases are 
of most interest for applications. Other values of the extension parameter correspond to 
additional contact interactions [30 a , and some of the values are of physical interest as well. 
To find a closed form of Green functions for the arbitrary value of the extension parameter 
is a more complicated task. The spectra of the corresponding extensions in the critical 
subspace are no longer periodic for such a situation that requires to apply more exquisite 
calculation methods. We suppose to consider this issue in our future publications. 

In addition, the nonrelativistic Green functions are constructed. The latter Green 
functions are represented for all possible types of 2+1 dimensional nonrelativistic particles. 
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6 Appendix 

For the sake of completeness we consider here the Green functions for the scalar particle. 
They are defined by Eqs. (J3J), (JH), ©, and ©, where S T (x,x') read 

S T (X, X) = ±i ^ ±0n (s) ±<fr* n {%') , 

and ±(p n (x) form a complete set of orthonormalized solutions of the Klein-Gordon equa- 
tion. Here we consider the natural extension of the Klein-Gordon operator for which 
solutions in 2 + 1 dimensions and the related spectrum read |3] , 

1 -i±ex\ jl J(l-loh 



± m , ; (x) = -7^ £ e- l±£X J^ l -^I m+H , lm (p) 



±e = ±VM 2 + u, oj = 7 [1 + 2m + \l + n\ + £ (I + fj,)] , 
I = 0,±1,±2, m = 0,l,2,.... 

Using Eqs. (j2*5Jl . P^jl. (|2*5jl . and (|50|l . we calculate the causal and anticausal propagators. 
They have the form 

S c (x, x') = / f sc (x, x', s) ds, S~ c (x, x') = / f sc (x, x', s) ds , 
Jo J-o 

r (x, x', s) = J2 fi C ^ *) . /r s) = A (s) e ^ e -W e -^j M > 



/ 

-ifieBs 



e 2 ^(z) 



r c (x,x', S ) = A( S )e- 
+F (z, Ay? - eBs, fi)+Y (z, -A(p + eBs, -//)] , (58) 

where A (s) is given in (J3T|) . and Y(z,r],fi) in (jHSJ. The expression (|HH)l can be 

generalized for the (-D + l)-dimensional case, where D is the number of spacial dimensions, 
with the substitution A (s) in (J5EJ) by v4 (Z)) (s), 

r . D W _izr\(0-2)/2 

^ oo = ^ oo ex p { h E ( A ^) 2 } ( ^ J > D ^ 3 • 
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